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Abstract. We examine the ability for the Large Area Telescope (LAT) to constrain 
Minimal Supersymmetric Standard Model (MSSM) dark matter through a combined 
analysis of Milky Way dwarf spheroidal galaxies. We examine the Lightest Supersym- 
metric Particles (LSPs) for a set of ~ 71k experimentally valid supersymmetric models 
derived from the phenomenological-MSSM (pMSSM). We find that none of these mod- 
els can be excluded at 95% confidence by the current analysis; nevertheless, many lie 
within the predicted reach of future LAT analyses. With two years of data, we find 
that the LAT is currently most sensitive to light LSPs (mLSP < 50 GeV) annihilating 
into r-pairs and heavier LSPs annihilating into bb. Additionally, we find that future 
LAT analyses will be able to probe some LSPs that form a sub-dominant component 
of dark matter. We directly compare the LAT results to direct detection experiments 
and show the complementarity of these search methods. 

Keywords: dark matter experiment, dark matter theory, dwarf galaxies, supersym- 
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1 Introduction 

Astrophysical evidence suggesting that non-baryonic dark matter (DM) comprises 
nearly 25% of the energy density of the Universe is one of the most compelling ar- 
guments for particle physics beyond the Standard Model (SM) [1]. At present, ex- 
perimental tests of this DM component are almost exclusively limited to gravitational 
interactions, and few constraints exist on the character of DM. Axions, dark pho- 
tons, sterile neutrinos and even more exotic theoretical constructs are all plausible DM 
candidates [2], though models containing a new neutral and stable weakly-interacting 
massive particle (WIMP) of mass ~ 100 GeV are by far the most studied. WlMPs are 
a favorable candidate because their mass and couplings to the SM can naturally give 
a cosmological relic density in agreement with the experimentally measured value [3]. 
Additionally, WlMPs point to new physics at the weak scale (~ 100 GeV — 1 TeV), a 
scale that has been the focus of much theoretical work to explain the stability of the 
Higgs potential and the origin of electroweak symmetry breaking. 

In the past several decades, supersymmetry (SUSY) has been the most widely- 
studied, and arguably the best-motivated, theoretical framework for physics beyond 
the SM [4-10]. In the most attractive SUSY models, an extra matter parity ("R- 
parity" ) symmetry is used to simultaneously explain the stability of the proton and of 
the Lightest Supersymmetric Particle (LSP). In viable SUSY models the LSP is often 
the lightest neutralino (Xi); which is one of the most widely studied examples of WIMP 
DM. Generic predictions of SUSY are difficult to obtain, since the minimal consistent 
SUSY extension of the SM, the Minimal Supersymmetric Standard Model (MSSM) 
introduces more than 100 free parameters. A typical strategy for overcoming this diffi- 
culty is to highly constrain this set of parameters by employing aesthetic assumptions 
about the physical origin of SUSY at a very high energy (i.e., mSUGRA [11, 12]). 
In contrast, here we study a broader and more comprehensive subset of the MSSM, 
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the phenomenological-MSSM (pMSSM) [13]. The pMSSM is derived from the MSSM 
using experimental data to ehminate parameters that are free in principle, but highly 
constrained by observations (e.g., sources of flavor violation in the new physics flavor 
sector). Thus, the pMSSM provides a compromise between the need to remain flexible 
and somewhat agnostic in assumptions about yet-undiscovered physics and the need 
to categorize the range of predictions made by well-motivated models. The LSPs of 
the pMSSM are viable candidates to comprise some or all of DM, and they may be 
probed through a variety of experimental approaches. 

The possibility of DM-SM interactions having weak-force strength allows an ex- 
citing opportunity to detect and characterize the nature of DM via a combination of 
experimental efforts. For example, weak-strength interactions might lend themselves 
to study at the LHC, where DM particles could be produced and studied indirectly 
through missing energy signatures. Additionally, the DM halo permeating our galaxy 
could be detected directly through scattering interactions between DM particles and 
nuclei in detectors on Earth. Finally, indirect detection of DM is possible through 
the astrophysical observation of anomalous energetic SM particles resulting from DM 
particle annihilation (or decay). 

One of the most sensitive instruments for the indirect detection of DM is the 
Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi). 
Gamma rays from the final state of DM annihilation (or decay) would be produced 
preferentially in regions of high DM density and may be detectable by the LAT. Dwarf 
spheroidal satellite galaxies (dSphs) of the Milky Way are promising targets for the 
detection of such a signal. These dSphs are DM-dominated and lack active astrophys- 
ical production of 7-rays [14, 15], a troublesome background in many other searches 
for DM annihilation. The LAT Collaboration recently presented results constraining 
the annihilation cross section for a small set of prototypical DM models from a joint 
likelihood analysis of 10 dwarf spheroidal galaxies [16]. In the present paper, we extend 
this analysis to an investigation of ~ 71k pMSSM models previously discussed in the 
literature [13]. 

We begin by briefly discussing the techniques employed to generate ~ 71k pMSSM 
models and the various constraints imposed in their selection. We next describe the 
combined likelihood procedure for setting upper limits on the annihilation cross section 
for each pMSSM DM model using LAT observations of ten Milky Way dSphs. We 
compare the LAT cross section limits to the actual cross section for each pMSSM 
model and study the SUSY model dependence of these results in detail. The main 
findings are: (i) that the LAT is currently most sensitive to light LSPs {m.^a < 50 GeV) 
annihilating primarily to r-pairs, (ii) that annihilations to r-pairs are actually harder 
to limit than annihilations to the other channels for relatively heavy [m^o > 50 GeV) 
LSPs, and {Hi) that, surprisingly, a significant fraction of the models that are near LAT 
sensitivity have LSPs that would form a sub-dominant component of the total DM halo. 
Additionally, we discuss the relationship between LSP eigenstate composition and the 
LAT sensitivity. Finally, we compare expectations for near-future LAT dSph searches 
and direct detection experiments. 
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2 Generation of the pMSSM Model Set 



The sensitivity of the LAT detector can be explored over a broad region of supersym- 
metric parameter space. We investigate the 7-ray production from ~ 71k points in 
the 19-dimensional pMSSM parameter space generated in previous work [13]. These 
points pass all of the constraints discussed in this section and are referred to as pMSSM 
models. 

The 19-dimensional parameter space of the pMSSM results from imposing the fol- 
lowing minimal set of assumptions on the general R-Parity conserving MSSM [17]: (z) 
the soft parameters are taken to be real, allowing no new CP-violating sources beyond 
those in the CKM matrix; (ii) Minimal Flavor Violation [18] is taken to be valid at 
the TeV scale; [iii) the first two generations of sfermions having the same quantum 
numbers are taken to be degenerate and to have negligible Yukawa couplings; and (iv) 
the LSP is taken to be the lightest neutralino and is assumed to be a stable thermal 
WIMP. No assumptions about the physics at high energy scales or SUSY-breaking 
mechanisms are employed. The first three conditions are applied to avoid issues as- 
sociated with constraints from the flavor sector. These assumptions reduce the SUSY 
parameter space to 19 free soft-breaking parameters that are given by the three gaug- 
ino masses, Mi o 3, ten sfermion masses rrin n ;i ;i f r z ^ , the three A-terms 
associated with the third generation (Abf ,-), and the usual Higgs sector parameters /i. 
Ma and tan/3. 

The set of models discussed in this paper was selected by numerical scans over 
the 19-dimensional parameter space of the pMSSM. This selection required a choice of 
parameter range intervals and scan priors. Issues involved in this selection have been 
described in detail previously [13, 19, 20]. Here, we simply note that two scans were 
performed: one employed a fiat prior beginning with lO'' points, and a second used a 
logarithmic prior employing 2 x 10® points. The relevant differences between these two 
scans are that (i) all SUSY mass parameters were restricted to be < 1 TeV for the 
fiat-prior case, while for the log- prior case this restriction was raised to < 3 TeV, and 
(a) the choice of the logarithmic prior generally leads to more compressed sparticle 
spectra than does the fiat-prior case. Note that the restriction on the upper limit for 
the mass parameters in both scans is chosen to ensure relatively large production cross 
sections at the LHC. In the present work, we focus primarily on results for pMSSM 
models in the fiat-prior set. 

After scanning the 19-dimensional parameter space, we subjected the resulting 
points to a set of theoretical and experimental constraints to select models that are valid 
for study. We briefly review these restrictions here.^ (i) Our theoretical constraints 
required that spectra must be tachyon free, color and charge breaking minima must be 
avoided, and a bounded Higgs potential must exist (leading to radiative electroweak 
symmetry breaking), (ii) We employed a number of constraints from the flavor sector 
and precision electroweak data arising from the measurements of {g — 2)^, b — sj, 
B — 7- rz/, Bs — fJ''^fJ'~, meson-anti-meson mixing, the invisible width of the Z and Ap. 
{iii) Restrictions resulting from numerous direct searches at LEP for both the SUSY 

^For full details, see ref. [13] 
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Figure 1. Logarithm of the LSP rehc density, Qlsp^^j for the flat-prior selected pMSSM 
models (grey bars) and log-prior selected pMSSM models (black bars). 

particles themselves, as well as the extended SUSY-Higgs sector, were imposed. Some 
of these searches needed to be re-evaluated in detail to remove SUSY model-dependent 
assumptions [13]. (iv) Null results from the set of Tevatron SUSY sparticle and Higgs 
searches were imposed. The most restrictive Tevatron data came from searches for 
stable charged particles [21] and from searches for an excess of multijet events with 
missing transverse energy [22]. We note that in the latter case, the search strategies 
were designed for kinematics expected in mSUGRA-inspired models, which required 
specialized simulations to apply these results in the context of the pMSSM [13]. While 
these bounds have been superseded by LHC analyses, we note that collider searches 
have been found to be largely uncorrelated with the most important DM observables 
(relic density, annihilation and scattering cross-sections) [23], and we expect the results 
presented here to accurately represent the prospects for indirect detection in the context 
of the pMSSM. (v) Finally, we have required that the LSP contribution to the dark 
matter relic density not exceed the upper bound determined by WMAP [3] and that 
the LSP scattering cross sections obey concurrent constraints from direct detection 
experiments. 

After imposing theoretical and experimental constraints, ~ 68.4k models from 
the fiat-prior sample and ~ 2.9k models from the log-prior sample remain for study. 
Note that the LSPs of the pMSSM models are not required to saturate the measured 
relic density. Thus, many models allow for multicomponent DM, a prosaic example 
being composed in part by a pMSSM neutralino and in part by an axion that solves 
the strong CP problem. The distributions of LSP relic densities found in our fiat- and 
log-prior model sets are shown in Figure 1. We find that a small subset of pMSSM 
models saturate the WMAP bound (with ^Lsp/i^ > 0.1). 

Gamma-ray energy spectra are calculated for all ~ 71k pMSSM models using 
the computational package DarkSUSY 5.0.5 [24]. DarkSUSY calculates the total 7- 
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ray yield from annihilation, as well as the rates into each of 27 final state channels. 
We omit contributions from the loop-suppressed monochromatic channels 77 and 'yZ^ 
which, although distinctive, are typically tiny in our model set. These 7-ray spectra 
are tested for consistency with the LAT 7-ray data. 

3 LAT 7-ray Limits on the pMSSM 

The Milky Way dSphs are a promising set of sources for the indirect detection of 
DM via 7 rays. Stellar velocity data from these galaxies suggest that they are very 
rich in DM, while observations at other wavelengths show no signs of astrophysical 
signals [14, 15]. The integrated signal flux at the LAT, 0^ (ph cm^^ s^^)) from pair 
annihilation in a DM distribution with density given by p{r) is 

The preceding J-factor represents the line-of-sight integral through integrated over 
a solid angle AQ, while the second factor is strictly dependent on the particle physics 
properties of the DM model and is integrated over the experimental energy threshold 
(discussed in more detail in section 4.2). Using Eq. (3.1), it is possible to combine the 
integrated J-factor for a dark matter distribution (calculated empirically for the dSph) 
with the particle physics characteristics of a DM model (from the pMSSM model set) 
to describe the predicted 7-ray annihilation signal. 

We follow the procedure of Ackermann et al. (2011; henceforth All) to constrain 
the 7-ray signal from ten dSphs with a joint likelihood analysis of the LAT data. Our 
data sample and event selection are identical to those described in All, taking photons 
in the energy range from 200 MeV < E < 100 GeV. In concordance with All, we use 
LAT ScienceTools^ version v9r20p0 and the P6_V3_DIFFUSE IRFs.^ J-factors and 
associated uncertainties for the ten dSphs are taken from Table 1 of All, where they 
were calculated using line-of-sight stellar velocities and the Jeans equation [16]. 

Our procedure for constraining (av) differs from that of All in that we model the 
7-ray emission from the dSphs with spectra generated from the ~ 71k pMSSM models 
rather than prototypical annihilation channels (i.e., bb, t~^t~ etc.). We calculate a joint 
likelihood for each pMSSM model by tying the pMSSM model parameters across the 
regions of interest (ROIs) surrounding the ten dSphs. Following All, we incorporated 
uncertainties in the J-factors of the dSphs as nuisance parameters in our likelihood 
maximization. Thus, our joint likelihood function is 

L{D I p^, {pk}) = l[L]:^^{Dk I Pn., Pk) 

(3.2) 

X I g-(logio(Jfc)-logio(Jfc))V2«T2 

ln(10)Jfcv^o-fc 

^http://ferini. gsfc.nasa.gov/ssc/data/analysis/software 

■^http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/overview. html 
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Here, k indexes the ROIs, L]^ denotes the standard LAT binned Poisson hkehhood for 
the analysis of a single ROI,^ represents the binned 7-ray data, {pm} represents the 
set of ROI-independent pMSSM model parameters, and {pk} are the ROI-dependent 
model parameters. Included in {pk} are both the flux normalizations of background 7- 
ray sources (diffuse and point-like) and the associated dSph J-factors and uncertainties. 
We find no significant 7-ray signal from any of the dSphs when analyzed individually 
or jointly for any of the pMSSM models. 

For each of the ~ 71k pMSSM models, we calculate the maximum annihilation 
cross section, (cif )ul, consistent with the null detection in the composite LAT data. 
To deal with nuisance parameters present in the joint likelihood, we constructed a 
profile likelihood Lp{D\{av)) by scanning in (av) and maximizing L(D | pm, {pi}) 
with respect to the other free parameters [25] . For each model, we proceeded to obtain 
a 95% one-sided confidence interval on the value of (av) by first maximizing log(Lp) 
with respect to {av) and then evaluating log(Lp) at increasing values of (av) until 
AlnLp = —2.71/2 [25]. This one-sided 95% confidence limit on (erf) serves as our 
value of (o"y)uL, which is compared to the true value annihilation cross section for each 
pMSSM model. 

4 Results 

In this section, we investigate the ability to constrain the predicted flux spectra of 
pMSSM models using LAT 7-ray flux measurements from ten dSphs. We then discuss 
the SUSY model dependence of our results and improvements that can be made in 
searches such as this one to enhance sensitivity to annihilations from certain classes of 
models. 

4.1 Model Constraint Distance 

We present limits placed on each of the ~ 71k pMSSM model spectra from a joint 
likelihood flt of ten Milky Way dSphs as described in Section 3. For each model, we 
compute the effective distance from constraint as the ratio 

V = {av)uL/{crv)T, (4.1) 

which compares the maximum annihilation cross section allowed by the data for a given 
spectral shape, (o"f)uL, to the actual cross section predicted for the pMSSM model 
giving rise to this spectral shape, {<jv)t- Here, the symbol (o"f )t represents the current 
annihilation cross section obtained assuming that the LSP follows a standard thermal 
cosmological evolution.^ It is important to note that, because we have employed the 
WMAP measurement of the DM relic density only as an upper bound when selecting 
pMSSM models, the appropriate normalization of the 7-ray spectra for each model 

^http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/Cicerone_Likclihood/ 
^One should not confuse the "T" in {(jv)t, here denoting thermal cosmology, for the thermal 

average over the current halo Boltzmann distribution, which is here denoted by the angle brackets 

(...). 
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involves a rescaling of the empirical estimate for the local DM energy density by the 
factor 

^=?7^' Mt^M^?', (4.2) 

i iWMAP 

where ^wmap^^ = 0.1143 [3]. The factors of R in {(Jv)t serve to obtain the cor- 
rect thermal number density of WIMPs n = p^o/m^o = {pQR)/m^o (where p^o and 
Po ~ 0.3GeVcm~^ are the energy density of the LSP and the empirically determined 
[26] total local DM energy density, respectively). LSPs in models with i? <C 1 are inter- 
preted as comprising one component of a multicomponent DM halo. For definiteness, 
we identify pMSSM models with LSPs that nearly saturate the WMAP relic density, 
^LSP^^ > 0.10 {R > 0.875), and distinguish this subset of models in the figures that 
follow. Deviations from the scaling in Eq. (4.2) would be appropriate for models where 
non-standard cosmology or other non-standard mechanisms [27-33] sever the relation- 
ship between the relic density and the present-day annihilation cross section obtained 
from thermal Boltzmann cosmological evolution. We assume that standard thermal 
evolution occurs and calculate results relative to this scenario. 

The ratio V (Eq. (4.1)) is calculated for each pMSSM model and is displayed 
in Figure 2. Models with V < 1 would be excluded at > 95% confidence, while 
models with V > 1 evade this limit by a factor of V. None of the pMSSM models 
are excluded at 95% confidence by this analysis; however, values of V reach ~ 1.5 for 
many of the models with LSP masses < 40 GeV, relic densities nearly saturating the 
WMAP measurement {R ~ 1) and annihilating predominantly to r-pairs (as is further 
discussed in section 4.2). 

Since the LAT is presently very close to constraining part of the pMSSM param- 
eter space, it is useful to estimate how constraints may improve over a 10 year mission 
lifetime. In the low-energy, background dominated regime, the LAT point source sen- 
sitivity increases as roughly the square-root of the integration time. However, in the 
high-energy, limited background regime (where many pMSSM models contribute), the 
LAT sensitivity increases more linearly with integration time. Thus, 10 years of data 
could provide a factor of y/b to 5 increase in sensitivity. Additionally, optical surveys 
such as Pan-STARRS and the Dark Energy Survey could provide a factor of 3 increase 
in the number of Milky Way dSphs corresponding to an increased constraining power 
of 73 to 3 [34]. On going improvements in LAT event reconstruction, a better under- 
standing of background contamination, and an increased energy range are all expected 
to provide additional increases in the LAT sensitivity. Thus, we find it plausible that 
the LAT constraints could improve by a factor of 10 compared to current constraints, 
and we choose to examine pMSSM models with D < 10 in detail. The interest in 
models with D < 10 is additionally motivated by predictions that DM substructure 
may increase the J-factors of dSphs by a factor of 2 to 10 [35, 36]. Such a boost would 
translate into a tightening of the current upper limits on (erf) by a comparable factor. 

In Figure 2 we observe, as expected, that the relic density of a given pMSSM 
model is important in determining the 7-ray signal strength. Despite this, we see that 
that the range of predictions in either subset of models is quite large, spanning many 



-7- 



10000. 
1000. 
•3 100.- 

o 

« 10.- 



2 4 6 8 10 12 -12 -10 -8 -6 -4 -2 2 

Logio[£)] Logi(,[<o-v>7- (10"^''cmV)] 

Figure 2. Distance from constraint, P (Eq. (4.1)), for all flat-prior selected pMSSM models 
(grey bars) and flat-prior models with i? ~ 1 (blue bars). We similarly histogram {av)T and 
note that the T) distribution is largely driven by the {(Jv)t distribution. 

orders of magnitude. We note the surprising fact that, among models with T) < 10, 
nearly 40% do not saturate the WMAP bound. 

When we compare our set of log-prior models to the flat-prior set, we find that 
the results are generally consistent with a somewhat larger tail towards large values of 
V. As noted in the next section, the shape of the V distribution is largely determined 
by the shape of the {(tv)t distribution, and one can verify that the difference in T> 
distributions between flat- and log-prior cases is echoed in their {av)T distributions. 
Models with such large V values are seen to be special cases, involving a finely tuned 
relationship between the DM mass and the mass of one or more of the SUSY Higgs 
states (such that m^o ~ rnhjl, so-called "Higgs funnel" models), that annihilate so 
efficiently in the early universe as to have minuscule relic density. 



4.2 SUSY Model Dependence 

In discussing the SUSY model dependence of these results, we first note (c.f.. Figure 
3) that most of the span in T> arises from the wide ranging values of the quantity 
((jf)T/2m-o. For a given value of (o"f)T/2m-o we observe that there is only about an 
order of magnitude span in D. To understand this finding we note that the total signal 
7-ray flux, as calculated in Eq. (3.1) (and repeated here for clarity), is 



10000. 
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Figure 3. Flat-prior pMSSM models represented in the D vs. {ov)t j'^vn^-o plane. Grey 
points represent generic models in this set while the subset of models with i? ~ 1 are high- 
lighted in blue. One observes that the wide range of T) values corresponds directly to the 
wide range ((Tt))r/2m?o values in our model set, and that at a given value of (ctu)t/2?tt,-o 

Xi Xi 

there is about an order of magnitude span in D values. It can also be noted that only ~ 60% 
of models with P < 10 have i2 ~ 1. 



We focus now on the particle physics-dependent piece, which is itself a product of the 
factor (cri;)T/2m|o and of the integral over the the total 7-ray continuum yield curve:^ 

dE^ ^ \ dE^ dE^ dE^ /■ ^ ' ' 

The V value expected from a given model is thus predominantly determined by the 
total thermally-averaged annihilation cross section, {(Jv)t, the relic density (via R) 
and the LSP mass (via (af )T/2m|o), and to a much lesser extent by the shape of the 
SUSY signal spectrum. 

We next discuss the SUSY model dependence of the spectral shape of the anni- 
hilation signal. It is useful to remove the large prefactor, (cxf )t/^-o, focusing on the 

Xi 

quantity: 

<Ps X {m\o/{(Tv)T) oc / ^^^^^7 > (4-5) 

where we integrate the total continuum yield curve over the experimental energy range 
from Emin = 200 MeV to E^ax = 100 GeV. In translating to 0^ x (m-o/ {(tv)t), we have 

Xi 

removed the explicit dependence on the total annihilation rate {av)T and most of the 



^ Here the sum is over annihilation final-state channels with terms describing hadronization yield 
("secondary" 7 rays), final-state radiation (FSR) and virtual internal bremsstrahlung (VIB). As pre- 
viously noted, a possible monochromatic 7 contribution is negligible here. We use the language of [37] 
in discriminating FSR and VIB, although it has been pointed out that such a distinction is somewhat 
artificial (or not even gauge invariant) [38]. 
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dependence on the LSP mass. What remains is a quantity that depends on the relative 
annihilation rates into SM final state channels, which is useful for comparing models 
with similar {av)T and LSP mass. A small residual LSP mass dependence remains 
in (ps X (m|o/(c'"f )t), via the relationship between the LSP mass and the limits of 
integration in Eq. (4.5). 

The supersymmetric origin of our 7-ray spectra is reflected in the distribution of 
annihilations into distinct SM final states (/). In SUSY, due to the Majorana nature 
of the annihilating Xi particles and the fact that x^s are non-relativistic in current 
DM halos, the annihilation rates for processes such as XiXi ~^ ff are proportional to 
{rrif /m^o)"^ , a fact that is often referred to as helicity suppression [39]. The ratio of 
rates into distinct channels (say // and /'/') is thus oc (mj/mj/)^ so that annihilation 
rates into channels with heavy final-state SM particles often dominate those with 
lighter final-state SM particles. These ratios are complicated by the fact that the 
rates XiXi ~^ f f are coherent sums of subprocesses with varying couplings, mixing 
angles, and exchanged particle masses. Nonetheless, in our pMSSM model set, we 
see that the dominant annihilation final states are often in accord with the helicity 
suppression intuition. For the vast majority of models in our set, annihilations are 
predominantly composed of a mixture of the bb, t~^t~, W~^W~, Z^Z^ and ti channels 
that are kinematically allowed (i.e., m^o > rrif). We observe a number of cases where 
the loop-level annihilation to gluons, XiXi ~^ 99) is dominant^ and a small number 
of cases where there are sizable (but still sub-dominant) contributions from the hA^, 
HA^, hZ^, HZ'^, W^H^ and (monochromatic) 77 or 7^° channels. 

The ability of the LAT to constrain a given pMSSM model can be decomposed 
into two planes (displayed in Figure 4). The 0ul (LAT fiux upper limit) vs. LSP 
mass plane demonstrates the ability of the LAT to constrain the spectral shapes of 
various annihilation channels, while the (ps x (^|o/ {o'v)t) vs. LSP mass plane describes 
the (scaled) signal fiux expected from each channel. We see that subsets of models 
with nearly pure annihilation final state channels cluster tightly in both the 0ul vs. 
LSP mass and cps x (m-o/ {(tv)t) vs. LSP mass planes and, with the exception of the 
dominantly channel models,^ extend down to their kinematic endpoints.^ 

^Loop-level annihilation to glouns can occur when the LSP co-annihilates with light flavored 
sfermions (e.g., m, d, e, etc.) in the early Universe. In current DM halos, annihilation by ex- 
change of the lightest sfermions is heavily helicity suppressed mj/rn^o) so that the loop- level 

process XiXi ~^ 99 may become the most efficient annihilation channel today. 

®While the subset of models with {{av) z z / {(^v)) > 0.42 do not extend down to m^^o w mz, the 
subset of models with {{av) z z / {crv}) > 0.20 do extend to this kinematic endpoint. The highest purity 
for XiXi ^ Z'^ Z'^ annihilations in our model set is {{<Jv) z z / {c'v)) ~ 0.45, such models also annihilate 
to X?X? W+W- with {{(Tv)ww I {crv)) w 0.55. 

^The set of 0(10) orange points near the top threshold, m-j^o « m^, are models that annihilate 
dominantly through the x5xi ~^ channel. They have very bino-like LSPs (supressing many other 
channels) and currently annihilate dominantly through stop exchange. These models satisfy the 
WMAP relic abundance constraint either by co-annihilation with a light stop or via the exchange of 
very light sfermions (i.e., channels that were more efficient at freeze-out). Since they are forced to 
annihilate dominantly to tt with m^o « nit, they are phase-space suppressed, {av)tt (1— /m^)^/^, 
and their fluxes are much lower than typical models that annihilate through this channel. 
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Figure 4. Flat-prior pMSSM models in the LAT flux upper limit vs. LSP mass (left panel) 
and scaled signal flux vs. LSP mass (right panel) planes (taking Jq = 10^^ GeV^ cm~^ as 
a reference J-factor). The full flat-prior model set is displayed as grey points and models 
whose annihilations occur predominantly through a given final state channel are overlaid in 
other colors. Models with [av) / {av) > 0.95 (green), with {av)^^^/ (av) > 0.93 (red), with 
{av)\y+\Y- / {(jv) > 0.95 (blue), with (av) z<izo / {cv) > 0.42 (magenta) and with {av)tt/{av) > 
0.85 (orange) are shown. Purities are chosen to obtain model subsets of similar size. 



The left panel of Figure 4 confirms that the LAT search places tighter constraints 
(0ul) on harder 7-ray spectra, which is to be expected as astrophysical backgrounds 
fall rapidly with energy. In particular, we see that LAT constraints on spectra from 
nearly pure XiXi ~^ t^t^ annihilations are tighter, by about an order of magnitude, 
than those placed on nearly pure XiXi ~^ bb annihilations. One expects the shapes 
and the relationship between these two "curves" to change significantly as LAT data 
taking continues. The LAT sensitivity to hard spectra, which contribute a significant 
number of photons in the background- free regime (> 10 GeV), is expected to increase 
more quickly than the sensitivity to softer spectra, which contribute in the background- 
dominated regime. Thus, the gap between limits on XiXi ~^ t^t~ and XiXi ~^ bb in 
the left panel of Figure 4 is expected to widen with increased LAT data taking. 

We note from the right panel of Figure 4 that nearly pure annihilations to r-pairs 
yield about an order of magnitude fewer signal photons than models with nearly pure 
annihilations to bb, at the same LSP mass and {(Jv)t- This finding was discussed at 
length in ref. [19], where it was noted that, although the 7-ray spectra resulting from 
r-like annihilations are harder (due to a large contribution from prompt vr'' decay), they 
are also shallower at low energies. This is demonstrated in Figure 5, where we display 
spectra (as calculated by DarkSUSY 5.0.4) for models that annihilate nearly purely 
into the bb and r-pair final-state channels. The curves displayed here correspond to the 
integrand in Eq. (4.5) and, when plotted in terms of the variable x = E.y/m^Q, have 
a nearly universal shape (there is significant SUSY model dependence in the signal 
spectra for ^ "^x?' internal bremsstrahlung [37]). Comparing the bb and 

r-pair cases it is clear that, at a given LSP mass and {crv)T, the softer bb spectra will 
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produce a much larger integrated flux of signal 7 rays than the r-pair spectra. We also 
observe that the vast majority of signal 7 rays have energy ^ ^x?' I'egcirdless of 
the particular final-state channel, so that the integrated signal flux depends sensitively 
on the lower limit of integration Xmm = Emin/fT^Xx- Heavier LSPs allow integration to 
lower X values and thus a wider gap between the total flux from hh and r-pair spectra. 
This is reflected in the relative slope of the red and green "lines" in the right panel of 
Figure 4. 

The two panels of Figure 4 describe two important factors in determining V = 
)uL/(o"f )r = 4>vl/4>s- (0 the ability to tightly constrain the flux from a given 
spectral shape (0ul) and (m) the flux that can be expected from that particular spectral 
shape {4>s)- These pieces of information must be combined to determine the ease with 
which the LAT can constrain final state channels (i.e., the magnitude of V), while 
providing a fair comparison between models annihilating to different final states with 
similar LSP mass and {(Jv)t- This is shown in the left panels of Figure 6. 

In the left panels of Figure 6 we display points for our models in the V x {av)T / mlg 
vs. LSP mass plane, with models that annihilate largely into single final state chan- 
nels colored as in Figure 4. Using V x {av)Tl'm?-Q is equivalent to using 0ul/(0s x 

Xi 

{mlo/ {o'v)t)), a ratio of the values found in either panel of Figure 4. We see that 
annihilations to a given final state are organized nicely in this plane. For comparison, 
in the right panels of Figure 6 we display points for each of our flat-prior models in 
the V vs. LSP mass plane using the same color scheme, where the colored points are 
seen to be highly mixed. Figure 6 allows us to predict whether it is easier to constrain 
annihilations to r-pairs or to bb (for example), at a given LSP mass, {(tv)t, and with 
experimental thresholds Emin and Emax- From the left panels we observe that annihila- 
tions to r-pairs are surprisingly more difficult to constrain than annihilations to bb for 
LSP masses m^o ^ 50 GeV, while the opposite is true for lighter LSPs. This crossover 
at m^o ^ 50 GeV is the point at which there is a balance between the relative ease of 
constraining the r-like spectral shape and the relatively low number of signal 7-rays 
produced in these annihilations, relative to the hb case. With longer observations, we 
expect this crossover to move toward higher energies as the limits on harder spectra 
are expected to tighten more quickly than limits on softer spectra. 

We emphasize that the grey points in the left panels of Figure 6 are as important as 
the colored bands. These regions represent the combination of annihilation final-state 
channels realized in the pMSSM model set and provide an estimate for the deviation 
from special cases of pure annihilations. For example, we observe that there is a special 
group of ~ 50 models with masses m^o 100 GeV that are constrained much more 
tightly than even the models annihilating purely to bb, i.e., the points falling below the 
red "line" in the top left panel of Figure 6. These are models (c.f.. Figure 5) that have 
a significant enhancement of the 7-ray spectrum near the endpoint x = E^/rriy^o f» 1. 
This enhancement is due to internal bremsstrahlung [40] [37], wherein an annihilation 
event produces a 7 ray directly from the hard process, i.e., XiXi ~^ ffl^ i^i addition 
to radiation off of the final state particles created from the process XiXi ~^ ff- The 
cross section for the 2 — )■ 3 process is naively much smaller than that for the 2 — 2 
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Figure 5. We display spectral shapes for pMSSM models that annihilate into bb (red) and r- 
pairs (green) final states with purities as in Figure 4. Spectra have been scaled to remove the 
((Tf)T/m|o prefactor and plotted in terms of the variable x = E^/m^o in order to emphasize 
the universality of individual final state spectra. 

process as the 2 — )• 3 case is suppressed by a fine structure constant and by three- 
body phase space. However, the leading terms in the 2 — )■ 2 process are helicity 
suppressed (~ (m//m^o)^) relative to the naive expectation, so that the 2 — )■ 3 process, 
XiXi ~^ ffl^ can be competitive with (or even dominate) the process XiXi ~^ //• This 
is especially true for annihilations to light final state particles (m^/m^o) ^ 1 through 
light superpartner mediators (m^/m^o) 1. Here, we find that many models in this 
group have SUSY mass spectra with very light sleptons, mg^ ~ m^o so that the would- 
be dominant annihilation channel XiXi ~^ e~^e~ is heavily suppressed and the 2 — )■ 3 
channel XiXi ~^ e'^e~'y contributes greatly to the overall 7-ray spectrum. 

4.3 LSP Eigenstate Composition 

The mapping of colored points from the left panels to the right panels of Figure 6, 
i.e., the combinations of ttl^o, (crv), R, and annihilation final-state distributions, is a 
complicated function of many of the couplings and masses that describe an arbitrary 
SUSY model. As such, it is difficult to robustly predict V for a given model by 
using only information about the interaction eigenstate composition of the lightest 
neutralino (our LSP DM candidate). While LSP eigenstate composition certainly has 
a large impact on the resulting 7-ray annihilation signal, it is difficult to disentangle 
the entire story from this piece of information alone. 

The impact of LSP eigenstate composition on the 7-ray signal is shown in Fig- 
ures 7-8. We display points similarly as in Figure 6 except that they are now highlighted 
in color according to whether their LSPs are dominantly bino, wino, higgsino or mixed 
(with purities as denoted in the caption). These figures show that, while the result- 
ing correlations are not as tight for eigenstate composition (Figures 7-8) as those for 
annihilation final-state (Figure 6), trends still exist. 
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Figure 6. We display points representing pMSSM models in the {av)ui,/m'io vs. LSP mass 

Xi 

(left panels) and T> vs. LSP mass (right panels) planes. The full flat-prior model set is 
displayed as grey points and models whose annihilations occur predominantly through a 
given final state channel are overlaid in other colors, as denoted in the figure. In the left 
panels, one can see that removing the dependence on total annihilation rate and LSP mass 
(scaling V x {av)T/m'io = (c^^)ul/'"^-o) allows for tight localization of models with similar 

Xi Xi 

annihilation spectra, whereas it is comparatively difficult to predict where models fall in the 
P vs. LSP mass plane without such scaling. The upper panels display these relations for all 
pMSSM models while the lower panels zoom in on those models that are closest to constraint 
{V < 10). 



Bino-like LSPs annihilate essentially only through channels XiXi ~^ ff via t- 
channel exchanges of sfermion partners /. Due to helicity suppression, bino-like models 
typically annihilate to some combination of t~^t~ , bb and tt (compare the red points on 
Figure 7 with the red/green/orange points in the top panels of Figure 6). The resulting 
mixture of rates into various final states depends on the pattern of sfermion masses, 
which are scanned over in our model generation procedure. The mapping in Figure 7 
from left- to right-panel is thus quite sensitive to scanned parameters, resulting in a 
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wide variety of predictions for 7-ray annihilation signals from bino-like LSPs.^*^ 

Wino-like LSPs annihilate dominantly to XiXi ~^ W^W~ by exchange of a 
chargino Xi that is almost degenerate in mass with the LSP.^^ These models are 
seen to correlate well with the case of annihilations purely to W^W~ (compare the 
blue points on Figure 7 with the blue points in the top panels of Figure 6). 

Higgsino-like LSPs couple efficiently to MSSM Higgs particles and to the 
annihilating through a multitude of channels regardless of the sfermion masses. When 
m^o > these models annihilate efficiently through XiXi ~^ Z^Z^ and through 
XiXi ~^ W^W~ via exchange of the second neutralino, ^'^^ the chargino Xi , 
which are nearly degenerate in mass with the LSP. Essentially all pMSSM models with 
higgsino-like LSPs annihilate dominantly through a mixture of the W^W~ and Z^Z^ 
final-state channels. The ratio of rates into the W^W~ and Z^Z^ final-state channels 
is primarily controlled by phase space (and asymptotes to {(tv) zz I {crv)ww ~ 0.75 for 
m^o mvF,niz).^^ The fact that these models annihilate to a mixture of W^W~ and 
Z'^Z^ final-state channels is refiected in the distribution of green points in Figure 8. 

In our model set, models with mixed LSP compositions are seen to essentially 
always have a substantial higgsino fraction and are generally described as a bino- 
higgsino mixture with a sub-dominant wino component. Given this, the mixed LSP 
distributions on Figure 8 are unsurprising. 

We observe that bino and mixed LSPs are over-represented in the set of models 
with "D < 10, as compared to their occurrence in the full model set. This is to be 
expected, since many such models are seen to achieve a large LSP relic density i? ~ 1 
while also maintaining a large annihilation cross section. In contrast, we find that 
none of the models in the P < 10 have LSPs that are classified as nearly pure winos or 
higgsinos. This is also in accord with our expectations as such purely wino or higgsino 
LSPs annihilate too efficiently in the early Universe, leaving i? ^ 1, at least for LSP 
masses below ~ 2 — 3 TeV. 

4.4 Comparison to Direct Detection 

It is expected that LHC searches, direct detection experiments, and indirect-detection 
experiments will provide highly complementary information about the nature of DM. 
The set of pMSSM models discussed in this work has already been studied in the 
context of LHC searches [20] [43], direct detection experiments [44] [23] and indirect 

^'^For a large number of models in our set the calculation of {<jv)t is also complicated by co- 
annihilations [41], wherein the relic density is set by an effective annihilation rate {av)eff in the early 
Universe, which is very different from the (av) governing annihilations in current DM halos. 

^^For models with very purely wino LSPs, the lightest chargino is highly degenerate in mass with 
the LSP. Null searches for new charged stable particles performed by CDF [42] and DO [21] have 
been applied in our model generation procedure [13], excluding models with stable charginos lighter 
than ^ 206 GeV. This results in a lower limit on the blue (wino-like) points (defined as having 
|Zi2|2 > 0.99) in Figure 7. 

^^Higgsino models with m^o < mw^'n^z typically annihilate mostly into r-pair or bb final states 
although we observe a number of models where the 77 and 7^° final states (proceeding through loops 
involving relatively light charginos) can become the dominant annihilation channels when the other 
relevant new particle masses are very heavy. 
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Figure 7. We display points in figures similar to those in Figure 6. Here the full flat-prior 
model set is displayed as grey points and models categorized according to LSP eigenstate 
composition are overlaid in other colors. By convention our LSPs are described in terms of 
their neutralino mass matrix entries as: Xi = ^iiB + Zi2W'^ + Zi^Hf + Z14H2. Bino models 
are defined as having > 0.99 and are displayed here in red. Wino models are defined 

as having |2^i2p > 0.99 and are displayed here in blue. 



detection searches with cosmic-ray electrons and positrons [19] and neutrinos [23]. We 
note that, essentially by construction, LHC searches are expected to rapidly exclude 
(or discover) most of the models in this set. Thus, a comparison of LAT and LHC 
results would be relatively unenlightening, since the most constraining LHC searches 
are typically only indirectly related to the detailed nature of the LSP. Therefore, we 
focus on a comparison between the prospects for future LAT dSph analyses and the 
limits on spin-independent and spin- dependent scattering cross sections. 

Since scattering signals are proportional to a single factor of the local dark matter 
density, direct detection experiments generally have an easier time setting limits on 
LSPs with low relic density (scaling like ~ i? as opposed to ~ for indirect detection). 
In Figure 9, we display the set of pMSSM models in the spin-independent ("SI," 
left panel) and spin-dependent ("SD," right panel) scattering cross section vs. LSP 
mass planes, highlighting the models within reach of future LAT dSph analyses (i.e.. 
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Figure 8. We display points in figures similar to those in Figure 6. Here the full flat-prior 
model set is displayed as grey points and models categorized according to LSP eigenstate 
composition are overlaid in other colors. By convention our LSPs are described in terms of 
their neutralino mass matrix entries as: Xi = ZnB + Z\2W'^ + Zi^H^ + Zi4i?2- Higgsino 
models are defined as having (|Zi3p + |Zi4p) > 0.99 and are displayed here in green. Mixed 
models are defined as having |-^i2p and (l^isp + |-Z^i4p) all < 0.7 and are displayed 

here in magenta. 

models with D < 10). We note that spin-independent scattering bounds have become 
significantly more constraining since the era when the pMSSM models were generated. 
The current best bound has been set by the XENONIOO Collaboration [45] and is 
depicted by the black curve on the left panel of Figure 9. An uncertainty of about a 
factor of four applies to this curve, due to uncertainty in the determination of matrix 
elements for nuclear scattering. 

We observe that there are many models that are expected to be discovered or ex- 
cluded by both direct and indirect detection experiments. This is a fortunate scenario, 
potentially allowing for relationships between LSP mass, (erf), R, annihilation final 
state channels and even details about heavier SUSY particles to be inferred. Addition- 
ally, we observe that there exist a number of models that will only be accessible to the 
LAT. These are models whose LSPs are dominantly bino and whose particle spectrum 
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Figure 9. Comparison of LAT dwarf search and direct detection limits. We display all 
points in the flat-prior model set in grey and points having I? < 10 in red. In the left 
panel the black curve depicts the current best SI scattering limit set by XENONIOO [45]. 
Near-future projected SI limits from LUX [47], SuperCDMS [48], COUPP 60kg and COUPP 
500kg [49] are displayed as, magenta-dashed, green-dashed, blue-dashed and blue-dotted 
lines, respectively. In the right panel Current SD scattering limits from the AMANDA 
[50] and IceCube-22 [51] collaborations are displayed as brown and orange lines, respectively 
(with the assumption of soft or hard channel annihilations represented by dash-dotted or solid 
lines, respectively). Near- future projected experimental limits from the COUPP [52] [49] 4kg, 
60kg and 500kg searches in blue- solid, dashed and dotted lines, respectively. The projected 
IceCube/DeepCore limit estimated in [53] is displayed as a magenta-dashed line (a more 
accurate IceCube/DeepCore analysis is presented in [23]). 

is somewhat hierarchical, including the light bino and one or more light sleptons. Such 
a scenario is essentially invisible in both direct detection experiments and at the LHC, 
due to a lack of accessible colored production channels. Generation of a new pMSSM 
model set that reflects progress in direct detection limits and the early running of the 
LHC is currently underway [46]. 



5 Conclusions 

We investigated the ability of the LAT instrument on board the Fermi observatory 
to detect supersymmetric dark matter. LAT observations of ten Milky Way dwarf 
spheroidal satellite galaxies were combined to constrain the 7-ray signal expected from 
a large set of phenomenologically-viable SUSY models. The LAT analysis sets some 
of the tightest constraints on annihilation cross section yet available from indirect 
detection experiments; however, it falls slightly short of ruling out any of the SUSY 
models in our set. We note that many models are quite close to being excluded 
(or discovered) by the LAT. We expect that future LAT observations and improved 
instrument performance will put significant pressure on this region of parameter space. 
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We found that the majority of models in our set have a SUSY-dependence in 
accord with intuition developed from highly constrained scans of the MSSM. However, 
many models have unique annihilation channels and rates that are only observed in the 
broader context of the pMSSM. We have investigated the relative abihty of the LAT to 
constrain annihilations into various final-state channels, noting that the LAT analysis 
is most sensitive to light LSPs {m^o < 50 GeV) annihilating dominantly to r-pairs and 
heavier LSPs annihilating dominantly to bb. Such behavior reflects a trade-off between 
the relative ease of constraining the spectral shape of annihilations to r-pairs with the 
relatively low number of 7 rays that are produced as a result of these annihilations. 
The sensitivity crossover point is expected to move toward higher energies as the LAT 
continues to take data. 

We have compared future expectations of the LAT dwarf search with those for 
direct detection experiments, finding examples of models that are accessible to combi- 
nations of the two experiment classes. Although the LAT search seems to be the more 
challenging method of discovering SUSY DM, we emphasize the unique character of 
such searches. Indirect searches are invaluable due to their sensitivity to DM signals 
regardless of the SM states that the DM couples to most strongly (rather than requir- 
ing a strong coupling to quarks and gluons as in direct detection experiments and, to 
a large extent, colliders). Additionally, the indirect detection of DM would provide 
information about the cosmological DM abundance, rather than having to infer DM 
properties from neutral detector-stabic particles produced at colliders. As the LAT 
mission continues, we expect it to extend the sensitivity of indirect searches for DM 
into this very interesting parameter space. 
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